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The mechanism for the 3x3 distortion of Sn/Ge(lll) 
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We show that two distinct 3x3 ground states, one nonmagnetic, metallic, and distorted, the 
other magnetic, semimetallic (or insulating) and undistorted, compete in a-phase adsorbates on 
semiconductor (111) surfaces. In Sn/Ge(lll), LSDA/GGA calculations indicate, in agreement with 
experiment, that the distorted metallic ground state prevails. The reason for stability of this state is 
analysed, and is traced to a sort of bond density wave, specifically a modulation of the antibonding 
state filling between the adatom and a Ge-Ge bond directly underneath. 



I. INTRODUCTION 

The study of the so-called a-phases (1/3 monolayer 
adatom coverage) of tetravalent adsorbates on semicon- 
ductor surfaces has recently attracted considerable in- 
terest, due to the complex and diverse phenomenology 
displayed by systems that, at a first glance, look very 
similar, both from a structural and from an electronic 
point of view. On one side we have Pb/Ge(lll) 0-0] 
and Sn/Ge(lll) |4|-[Tl|, where a transition from v^x v3 
to 3 x 3 surface periodicity has been observed below ~ 200 
K. On the other side, SiC(OOOl) Q and K/Si(lll):B Jl]| 
retain a V3 x v3 periodicity at all temperatures, but 
are insulating, in contrast with simple electron count- 
ing rules. All the above systems are characterized, in 
the \/3 x y/3 phase, by a narrow and half-filled sur- 
face band arising from the dangling bond orbital of the 
adatom. Narrow metallic bands are highly unstable ei- 
ther against electron-electron instabilities (Mott transi- 
tion), or against genuine structural distortions aimed at 
lowering the electronic density of states at the Fermi 
energy @. SiC(OOOl) and, possibly, also K/Si(lll):B 
appear to belong to the former class, due to the large 
Coulomb repulsion within the dangling bond orbital of 
the Si adatom. In line with this expectation, SiC(OOOl) 
has been recently predicted to be a surface magnetic 
Mott insulator based on first-principles calculations jl6) . 
Sn/Ge(lll) and Pb/Ge(lll) belong instead to the class 
of surfaces where a large structural distortion, leading to 
a 3 x 3 periodicity, removes at least some of the original 
destabilizing metallic character. 

Although considerable progress has been made in the 
last years, a complete understanding of the physics un- 
derlying the appearence of so different phenomena (Mott 
transition or atomic distortion) at the surface of other- 
wise very similar systems, is presently missing. In partic- 
ular, we lack a microscopic understanding of the reasons 
that make a system decide for one or the other state. 
This is particularly relevant also in connection with the 
possibility that these states may be observed at the sur- 
faces of other, presently unexplored, ev-phases such as 



Sn/Si(lll), Pb/Si(lll). Recent work on Sn/Si(lll) 
fails to indicate any sign of transition down to 100 K. 

Here we focus on Sn/Ge(lll), as a prototype of the 
systems where a large atomic distortion takes place. We 
investigate the system with first-principles methods, and 
find that both an undistorted (i.e. structurally y/3 x V3) 
but magnetic state, and a 3 x 3 structurally distorted 
state, lower the energy of the originally metallic y3 x \/3 
surface. However, the energy gain is larger for the struc- 
turally distorted case, explaining the observed low tem- 
perature 3x3 structure. Moreover, by examining in de- 
tail the atomic and electronic structure of the 3x3 distor- 
tion, we are able to highlight the microscopic mechanisms 
that drive the transition. 



II. METHOD 

The Sn/Ge(lll) surface has been modelled in a re- 
peated supercell geometry where three bilayer slabs of 
Ge atoms are separated by equivalently thick vacuum re- 
gions. Sn adatoms are placed in the T4 position of the 
upper surfaces while dangling bonds of the lower surfaces 
are saturated by hydrogen atoms. We performed exten- 
sive electronic structure calculations for both the V3x y/3 
and the 3x3 surfaces, either in the local (spin) density 
approximation (LSDA) or including gradient corrections 
(GC) |L8| to the energy functional. Norm-conserving 
pseudopotential |l^| in the Kleinman-Bylander form [^0| , 
a plane- wave basis set with 12 Ry energy cutoff, and a 
15 x 15 k-points grid to sample the full surface Brillouin 
zone (SBZ) of the \/?> x \/3 phase were used. For the 
3x3 surface an equivalent SBZ sampling was employed. 
All systems were structurally relaxed, keeping the bot- 
tom Ge-layer and the saturating hydrogen atoms fixed, 
until the Hellmann-Feynman forces on all other atoms 
were reduces to less than ps 10~ 2 eV/a.u. 
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FIG. 1. Surface band structure of the Sn/Ge(lll), calculated in the GC approximation. Left panel: \/i x \/3 undistorted 
nonmagnetic surfaces. Middle panel: 3x3 magnetic surface (here the structural distortion is negligible). Solid line, spin up 
electrons; dotted lines, spin-down electrons. Right panel: 3x3 distorted structure (nonmagnetic). In all cases shaded regions 
represent projected bulk bands. 



III. RESULTS 

The GC band structure of the unreconstructed Sn/Ge 
(111) -\/3 x -\/3 surface are reported on the left panel of 
Fig. [l| The system is metallic with a single predomi- 
nat, partially occupied, surface band, originating from 
adatom dangling bonds, in the bulk projected energy 
gap. Due to the rather small bandwidth, w ~ 0.6 eV, the 
system presents an instability toward a spin- and charge- 
density wave of 3 x 3 periodicity. The corresponding 
ground-state band structure is shown in the middle panel 
of Fig. |l|. In agreement with experimental evidence, the 
magnetically ordered structure is still metallic, although 
with a reduced density of states at the Fermi energy. A 
similar result was obtained earlier for the Si/Si(lll) sur- 
face ]2l[ where, however, the exchange gap was larger 
and the resulting structure was insulating. The inclusion 
in the calculation of GC terms to the energy functional 
turns out to be necessary to stabilize the magnetic struc- 
ture. This is not surprising as we expect the underling 
physics to be rooted in the same aspects (small band 
width, large coulomb self-interaction) that produce the 
Mott insulating states in related systems ||T^ , |l6|| . 

In addition to this magnetic instability metallic 
Sn/Ge(lll) \/3 x \/3 surface is also unstable against a 
purely structural distortion, where vertical displacements 
of the adatoms with 3x3 periodicity, are accompained by 
a bond alternation of the underlying substrate. In Fig. 
[| the atomic structures of the unreconstructed (left) and 
reconstructed (right) surface are compared. One of the 3 
adatoms rises above the surface while the other two sink 
deeper in the substrate, the final calculated vertical dis- 
placement between the two adatom types being as large 



as 0.36 A. The energy gain of the system upon recon- 
struction is of 9 meV/adatom, to be compared with 1-2 
meV/adatom of the magnetic case. The reconstruction is 
robust with respect to the exchange-correlation scheme 
used, and LDA and GC give here essentially the same 
results, indicating that the physics involved is probably 
well described by a conventional band picture, in spite of 
the large value of U/W in the surface state band. This 
point will be addressed later. Our calculated distortion 
compares favorably with that extracted from very recent 
x-ray diffraction jll|] and also with a previous LDA cal- 
culation 

The reconstruction develops with no energy barrier: 
we checked that the energy decreases steadily when the 
adatom vertical offset is fixed at about 10 % of its final 
value, the substrate atoms being allowed to relax from 
the unreconstructed positions. The role of the substrate 
relaxation is important: if only the adatoms positions 
are optimized, then the energy gain disappears, and the 
unreconstructed surface is more stable. Infact, it can be 
seen from Fig. || that the reconstruction pattern involves 
changes in the (111) bonds-lengths between second and 
third Ge atomic layers. The reconstruction apparently 
occurs because the symmetry lowering allows some re- 
hybridization between adatoms, accompained by a bond 
alternation in the substrate that stabilizes the deforma- 
tion. 

FIG. 3. Contour plot of the mid-gap state at M point in 
the SBZ of %/3 x y/3 surface. Full lines enclose regions where 
the wavefunction is positive; dotted lines enclose negative re- 
gions. Note the antibonding character of the wavefunction 
between the adatom and the substrate. Atomic positions and 
bonds are also shown for reference. 
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FIG. 2. Atomic geometrical structure of the Sn/Ge(lll) surface for the undistorted (y/3 x left panel) and the distorted 
(3 x 3, right panel) structrure. Atomic distances are in A 



The mechanism can be better understood by examin- 
ing the nature of the Wannier function (WF) of the mid- 
gap state of the V3 x V3 surface (approximated here by 
the wavefunction at the M point, see Fig. ||). The WF 
originates from the adatom dangling bond but has impor- 
tant contributions from substrate states. In particular it 
can be seen as an antibonding combination of the adatom 
dangling bond and the bonding state located in the (111) 
Ge-Ge bond below it. When the system reconstructs one 
of the three adatoms becomes inequivalent filling its state 
completely. By doing so the corresponding (111) Ge-Ge 
bond is strenghtened (due to the bonding character of 
the WF in that region) while the Sn-Ge one is weakened 
(the WF is antibonding there); as a result the adatom 
rises above the surface. The WF's centered on the other 
two adatoms are partially depopulated and an opposite 
relaxation occurs. 

To better illustrate this mechanism let us consider a 
toy system where the relevant structural motif (a Sn 
adatom and the four Ge atoms underneath) is extracted 
from the surface and hydrogen atoms are added to satu- 
rate Ge dangling bonds. The neutral cluster (Fig. |], left 
panel) corresponds to the unreconstructed surface (Fig. 
||, left panel) and presents a semioccupied highest molec- 
ular orbital (HMO) with antibonding character between 
Sn and Ge. 

When the HMO is completely filled, making the cluster 
negatively charged, all Sn-Ge bonds weaken and the cor- 
responding distances increase. Not much else happens; 
in particular the already strong vertical Ge-H bond is 
not modified significantly. The effect is more dramatic in 
the opposite situation, when the HMO is emptied mak- 
ing the cluster positively charged (Fig. [|, right panel). 
As expected from its antibonding character, depopulat- 
ing the HMO strenghtens Sn-Ge bonds shortening their 
distances, while the lower Ge-H bond is essentially de- 



stroyed. 

Why does Sn/Ge(lll) energetically prefer to distort 
rather than become Mott-Hubbard insulating ? Our re- 
sults suggest that the 3x3 distortion is based on the 
strong antibonding interaction between the adatom dan- 
gling bond and the Ge-Ge bond directly underneath. The 
fact that the energy gain should come by an alternating 
hybridization/dehybridization of the surface band with 
a deeper Ge-Ge bond has two implications. The first is 
that the large value of U/W inside the surface band is 
not very relevant to this state, unlike the Mott-Hubbard 
state. The second is that it shows that band-Jahn- Teller, 
a terminology often used to describe it, is not a cor- 
rect characterization for this state. It is rather a bond- 
density-wave. In this case the energy gain does not come 
from gap-opening, that is only partial, but from the mod- 
ulation of the strength of the Ge-Ge bonds under the 
adatoms. Large adatoms, narrow semiconductor gaps 
and a deformable lattice favor that. These conditions are 
not met e.g. in SiC(lll), where moreover poor screening 
enhances the value of electron repulsion (U). 

In conclusion, we have studied the competition be- 
tween magnetic and distorted nonmagnetic ground states 
of Sn/Ge(lll). Dominance of the latter has been under- 
stood as due to a modulation of the antibonding part- 
nership between the adatom and the underlying Ge-Ge 
bond. 

Our calculations were performed on the CINECA 
Cray-T3E parallel machine in Bologna, using the paral- 
lel version of the PWSCF code. Access to the Cray ma- 
chine has been granted within the Iniziativa Trasversale 
Calcolo Parallelo of INFM and the initiative Progetti di 
ricerca di rilevante interesse nazionale of MURST. Spon- 
sorship from INFM/LOTUS and from COFIN97 is also 
acknowledged. 
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FIG. 4. Ground state geometric structure of a SnGe4H7 cluster in the neutral (left panel), negatively charged (middle panel) 
and positively charged (right panel) states. Atomic distances in A. 
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